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1. INTRODUCTION 

The sandwich structure is constructed from a thick core material bonded to the skin materials on its 
upper and lower sides. The skin material usually has a high stiffness characteristic; this may be obtained by 
using fiber-laminate composites or high-modulus metals. The core structure can accommodate high 
compressive and shear strengths. Combining the skin and core will create a sandwich structure with the 
ability of high flexural modulus due to its high moment of inertia (insertion of material in the structure). 

Sandwich structures are commonly used in aeronautics, especially in aircraft (Boeing 787 uses 
composite sandwich panels as part of the wing skin structure [1]. During the operation, aircraft can fly at both 
low velocities during landing and approach, and high velocities when taking off and in cruise conditions. 
With this condition, aircraft are prone to impacting objects during flights, which may cause damage and even 
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lead to catastrophic events. Added to its disadvantages such as high manufacturing cost and longer 
manufacture time, we have an incomplete understanding of sandwich structures under impact loading [2]. 

Core material with a negative poisson ratio (NPR) may contract when it is compressed and expand 
when it is tensioned [3]. The term “auxetic” was first introduced in 1991 by Kenneth E. Evans [3]. Under 
low-velocity impact [4]-[11], high-velocity impact [12]-[15], and blast loading [16]-[20], a core with NPR 
performs better in terms of energy absorption, damage localization, and densification. With such an “auxetic” 
characteristic, the possibility of increasing the performance of sandwich structures is promising. 

Several authors have studied the impact behavior of sandwich structure with an auxetic core. An 
investigation finds the impact behavior of a composite sandwich structure, where the core is made of 
aluminum and the plate is made of steel, with auxetic and honeycomb core geometries under blast loading 
using finite element software (ABAQUS/Explicit). They found that the auxetic core was better at localizing 
the stress distribution, where the deformation concentrates towards the impacts area, than the honeycomb 
core due to the auxetic behavior [16]. 

Others investigated the low-velocity impact behavior of carbon-fiber-reinforced plastic (CFRP) 
sandwich panels with various geometries of auxetic and non-auxetic cores. They found that the auxetic core 
increases the impact resistance of sandwich composite due to core densification that reduces the indentation 
depth. The auxetic core effectively improved impact damage resistance at a higher impact energy level of 
76 J, whereas foam core performs better at energy levels lower than 10 J [4]. 

Studies from researchers worldwide indicate that there is a need and possibilities for better- 
performing structures for impact mitigations. Research has been done by implementing these auxetic unit 
cells as their numerical test subjects. The behavior of a 3D re-entrant structure under various crushing 
velocities. In their study, it is concluded that the auxetic of the structure may have greater energy dissipation 
when it is given a crushing loading. The energy dissipated increases as the crushing velocity are increased. 
Densification gives greater strength and higher crushing resistance performance [21]. 

Others also investigate that auxetic foam as a filler in a square tube under quasi-static compressive 
loading. These studies compare with conventional positive poisson's ratio (PPR) foam as a filler on the 
square tube. The result shows that square tubes with auxetic foam filler have higher energy absorption 
capability and specific energy absorption (SEA) than conventional PPR foam-filled square tubes. These 
studies found that auxetic structures can mitigate impact loading by densifying their geometry and have 
better energy absorption abilities through their auxetic characteristics [22]. 

This paper addresses technical notes on what is important in a finite element procedure to simulate 
sandwich plates with composite faces and a core with NPR under low-velocity impact. We used 
ABAQUS/Explicit to develop the models based on [4]. The aim of this work is to develop a model of a 
sandwich structure with an auxetic core, conduct numerical experiments of sandwich structure with an 
auxetic core under the low-velocity impact, and compare the result from Usta [4] which uses LS DYNA with 
this paper’s findings. The model is strictly limited to the development of a sandwich model under dynamic 
(impact) loading using ABAQUS/Explicit software. The materials were CFRP for the faces and polylactic 
acid (PLA) for the core. We strictly adopted the sandwich dimension from Usta [4]. 


2. RESEARCH METHOD 
2.1. Finite element method (FEM) modeling strategy 

The sandwich structure model was defined according to the geometry and material properties used 
by Usta et al. [4]. The dimension of the face is 150 x 100 mm using a 3D deformable shell with a thickness 
of 1.05 mm. The heterogeneous core resembles the same dimension in XY-plane and uses a 3D deformable 
shell with a total core thickness of 24 mm and 0.7 mm thickness for the auxetic cell struts. It is important to 
note that when creating the heterogeneous auxetic core, when the column height has non-zero decimals, it 
needs to be made directly as one (1) column and not one (1) unit cell and reproduced in the direction of the 
row and column. This strategy is to prevent errors when the geometry has meshed. The three-dimensional 
visualization of the extruded auxetic core is shown in Figure 1. 

The material used for the face is UD CFRP consisting of three plies with [0/45/90] orientation and 
utilizes the composite layup sub-module (see Tables 1 and 2). For the core, PLA polymer is used (see 
Figure 2 and Table 3). The impactor uses a 3D rigid shell with a hemispherical-end cylinder. The radius 
dimension of the hemispherical end is 10 mm and 35 mm in the length of the cylinder body. The impactor's 
mass is 22.5 kg, so the density and mass can adjust with the impactor's geometry as it is a loading media 
rather than the research interest. 

The constraint between the face and core was connected using "Tie Constraint" because we want to 
model the two components as if they are bonded together. Therefore an adhesive layer is not required to be 
modeled. However, if the users want to simulate the effect of an adhesive layer, then a different approach can 
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be made. For the master surface, the surface will be from the face. For the slave surface, the surface will be 
from the core. The surface of the face in which the constraint is defined is only one connected face-to-face 
with the heterogeneous auxetic core surfaces. The impactor is constrained as a rigid body with free 
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movements only in Z-axis translation (same direction as the loading direction). 


Figure 1. Re-entrant auxetic unit cell whole structure in a three-dimensional (3D) view 


Table 1. UD CFRP (elastic and damage properties for lamina [4] 


Parameters Symbols Value Units 
Elastic Properties (Engineering Constants) 
Density p 1.5x10-° — tonne/mm? 
Longitudinal stiffness E11 160 000 MPa 
Transversal stiffness E22 6500 MPa 
Out-of-plane stiffness E33 6500 MPa 
Poisson’s ratio Viz 0.12 - 
V13 0.12 - 
U23 0.40 - 
Shear modulus Giz 5680 MPa 
Gi 5680 MPa 
G33 2840 MPa 
Damage Properties 
Longitudinal tensile strength Xr 1020 MPa 
Longitudinal compressive strength Xç 650 MPa 
Transverse tensile strength Yr 70 MPa 
Transverse compressive strength Yc 100 MPa 
Longitudinal shear strength SL 80 MPa 
Transverse shear strength Sr 80 MPa 


Table 2. UD CFRP damage properties (fracture energy) for lamina [23] 


Parameters Symbols Value Units 
Longitudinal tensile fracture energy Gyr 48.4 N/mm 
Longitudinal compressive fracture energy Gye 60.3 N/mm 
Transverse tensile fracture energy Gyr 4.5 N/mm 
Transverse compressive fracture energy Gyc 8.5 N/mm 


25 


Stress (MPa) 
a 


5 


>O U 


PLA Stress-Strain 


0 0.01 0.02 0.03 0.04 0.05 0.06 
Strain (-) 


Figure 2. PLA stress-strain curve [4] 
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Table 3. PLA material properties [4] 


Parameters Symbols Value Units 
Polylactic Acid Material Properties 

Young’s Modulus E 1900 MPa 

Poisson’s Ratio v 0.36 - 

Density p 1.4x107-? — tonne/mm* 

Ductile Damage Properties 

Fracture Strain - 0.02398 - 

Stress Triaxiality - 0.333 - 

Strain Rate é 0 s! 

Displacement at Failure - 0.025892 Žž - 


The loads are defined by the impactor, a rigid body, and a reference point. The reference point is 
optional for load definition simplifaction, where it will couple all the impactor's structure, and the designated 
velocity will define only the reference point. The simulation is loaded under low-velocity impact with 2.6 m/s 
or 76 J. The boundary condition on the support implies the ASTM-7136 support configuration with a fixed 
support type (see Figures 3-5) for the boundary condition configuration). This boundary condition has been 
proven to avoid a ringing effect on the structure due to overconstrained boundary conditions, where sufficient 
support is used to hold the structure in place during impact. 


Figure 3. ASTM-7136 boundary condition configuration of simulation [4] 


(a) (b) 


Figure 4. Boundary condition location for (a) the upper face and (b) the lower face, as well as the partition of 
the part 


Figure 5. Final visualization of the boundary condition when applied to the model 
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The meshing uses an S4R element type where the mesh size is 1.25 mm with quad-structured 
elements. A partition of the face is recommended to connect the mesh between the face and the core. The 
impactor mesh is coarse as it has already been defined as a rigid body. 

The interactions use interaction property where the tangential behavior is defined with "Penalty" in 
the friction formulation section and a 0.2 friction coefficient. Also, "Hard" contact is chosen in the normal 
behavior section. The contact type is general contact (explicit) which includes "All" with self" so that when 
the structure deforms and contacts each other, they will notice neighboring elements and deforms with 
contact accordingly. 

The end time or the simulation is set depending on the units used. In our case, the time is 
0.015 seconds (or 15 milliseconds). The step increment can also be selected according to the user's objective. 
However, in our case, we will not change it and let it be default. The frequency interval of the data saving is 
also important to record how much data is and at what time the data is stored. The author uses space-time 
intervals evenly with the interval of 150 to get information for every 0.1 ms. The bigger the interval, the more 
data is recorded, and the more computer memory is needed. 

When using a composite material, to visualize the Hashin damage, damage initiation criteria 
(DMICRT) will need to be activated for the composite components so that the Hashin damage can be set as 
an output. The reason is that in ABAQUS/Explicit, the Hashin damage initiation criterion cannot be selected 
directly in field output. It can only be selected directly in the history output. 


2.2. Hashin damage initiation criteria 

We employed the Hashin damage criterion, a material damage initiation for fiber-reinforced 
composites. It is based on [8], [9], where the criteria consider four different damage initiation mechanisms, 
which are fiber tension, fiber compression, matrix tension, and matrix compression. In (1)-(4) show the 
equation for each damage criterion. Generally, when an initiation criterion has reached a value of 1.0 or 
higher, the damage initiation has been met, and the damage process will start on the structure. 


Fiber tensile 


(Gt) +0(%2) o) 
Fiber compression 

Ge) a 
Matrix tensile 

HE © 


Matrix compression 


2 ON 2 
a) F (S 1| yc T sh (4) 
Where, 011, 022 is the normal stress in direction-1 and -2. T43 is the shear stress in plane-1 and direction-2. 
XT, XC is the tensile and compressive strength in the fiber direction. Y7, Y° is the tensile and compressive 


strength in the matrix direction. S$“, ST is the longitudinal and transversal shear strength. a is the coefficient 
for the contribution of shear stress to fiber tensile initiation criteria. 


2.3. Explicit time integration 

The finite element method is used to solve differential equations, especially to find components such 
as stresses on a structure numerically. In specific, ABAQUS/Explicit uses explicit methods and a central 
difference method scheme. An explicit method is used as the model is dynamic loading, nonlinearity, and 
high discontinuity in the solution such as impact or failure simulation, especially in a three-dimensional 
model. In (5) is the explicit method (6) is the dynamic equilibrium (7) is the central-difference method 
equation. 


Y(t + At) = F(Y(t)) (5) 
Where, Y(t) is the current state and Y(t + At) is the future or time-step state, At is the time-step increment. 
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Where M represents the lumped mass matrix to calculate the accelerations of the nodals, P is the external 
load, and I is the internal load. 


l = Viri ia (7) 


ax. 2Ax 


3. RESULTS AND DISCUSSION 
3.1. Visualization of deformation 

The deformation visualization of the sandwich structure for whole and half structures is shown in 
Figure 6 and Figure 7. The damage progression of the sandwich structure is shown in Figure 8. The 
maximum deformation is located at the center of the sandwich structure, where the impact occurs. The 
phenomena of the structure do not visualize the correct behavior where the face made of composite is brittle 
and should damage upon impact. However, this phenomenon can differ when using different materials, such 
as metals. The composite material model needs to be re-evaluated. 


Figure 6. Sandwich structure deformation visualization for the whole structure 


Figure 7. Sandwich structure deformation visualization for half-structure 


3.2. Visualization of failure 

Hashin damage criteria visualization in fiber-tensile direction is shown in Figure 9. The progression 
of the damage is shown from 2.0 ms until 11.0 ms range with an increment of 1.0 ms. It can be seen further 
that some elements have reached the maximum value for the damage criteria (1.0) but have not yet been 
deleted (in default, damaged elements are deleted). The cause of these findings still needs to be investigated 
later. 

As the simulation time increases, the element also stretches larger than its original size, 
strengthening the leading hypothesis that there could be a possibility of missing or incorrect damage 
properties. Thus, the result does not show how composite panels under impact loading should be. In 
ABAQUS/Explicit, the Hashin damage model requires fracture energy properties which from Usta [4] is not 
available. Therefore, an effort is made to obtain the fracture energy properties. Some references refer to 
works by Koloor [23] for fracture energy. 
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The difference in software usage can be discussed as reference [4] uses LS DYNA as the software to 
investigate their work. However, in this paper, ABAQUS/Explicit is used as the main software. LS DYNA 
failure uses such failure criteria. Where a damage initiation point is reached, the element will instantly be 
deleted. This case is incomplete when conducted in the same manner in ABAQUS/Explicit, where damage 
evolution comes after the damage initiation part, which hence would not instantly delete the element after a 
damage initiation point is reached. The difficulty arises because damage evolution requires an important 
property: the material’s fracture energy. Incompatibility of all the material properties leads to incorrect and 
invalid simulation, thus misleading results in the damaged region. It is known that fracture energy is a 
material property usually obtained from experimental results. Therefore, an initial experimental test is needed 
for this material's properties to get a more precise result and data. 
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.00 
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t=3.0 ms 


Figure 8. Damage progression of the composite sandwich structure with thermoplastic auxetic core under 
low-velocity impact loading 
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Figure 9. Hashin damage—fiber tensile (DMICRT-HSNFTCRT) visualization of the composite sandwich panel 
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3.3. Plotting displacement, velocity, and acceleration 
The graph showing the displacement, velocity, and acceleration of the impactor is shown in 


Figures 10-12, respectively. This correlates with the basic principles of impact mechanics. Whereas the 


displacement increase, the velocity decreases, and the acceleration decreases, hold, and increase to reach 
Zero. 
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Figure 12. Acceleration graph of the impactor 


3.4. Plotting force-time and energy diagram 

The sandwich structure that the author's model is based on [4]. Therefore, the force-time and force 
displacement are compared between the author’s simulation and the reference data shown in Figure 13 and 
Figure 14, respectively. However, the energy diagram and the absorbed energy per component are not 
compared to the reference because the reference does not give information regarding the content. However, 
the author gives information regarding the total, kinetic, and absorbed energy shown in Figure 15. The 
absorbed energy per component is also shown in Figure 16. 

The result still needs re-evaluation as the behavior does not agree with the real-life behavior. It is 
discussed in section 3.2 regarding the difference between the simulation of [4] and the findings in this paper. 
Due to incomplete material property data in all the required sectors, certain obstacles are found, especially in 
the damaged region. Hence, the difference between this paper's results and reference [4] is acknowledged. 


Therefore, the objective of this visualization is to give a brief overview that the visualization can be obtained, 
despite the validity of the numerical value. 
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Figure 15. Energy diagram of the sandwich structure Figure 16. Energy absorbed per component of the 


4. 


sandwich structure 


CONCLUSION 
We provide a finite element procedure to simulate sandwich plates with composite faces and a core 


with an NPR under low-velocity impact. The model developed here could develop the simulation up to the 
elastic region. However, it could not fully resemble and represent the model from reference. Further 
improvements in the material modeling strategy, especially in the fracture modeling of the composite face 
with compatible material properties in all required sectors, especially damaged sections, are strictly 
necessary. 
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